Tumor Gompertzian growth by cellular energetic balance 
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A macroscopic model of the tumor Gompertzian growth is proposed. The new approach is based 
on the energetic balance among the different cell activities, described by methods of statistical 
mechanics and related to the growth inhibitor factors. The model is successfully applied to the 
multicellular tumor spheroid data. 



A microscopic model of tumor growth in vivo is still 
an open problem. It requires a detailed description of 
cellular interactions and a control on the large variety 
of in situ conditions related to the distributions of nu- 
trient, of oxygen, of growth inhibitors, of blood vessels 
and capillarity and to the mechanical effects due to tis- 
sue elasticity and heterogeneity. On the other hand, in 
spite of the previous large set of potential parameters, 
tumors have a peculiar growth pattern that is generally 
described by a Gompertzian curve Q , often considered as 
a pure phenomcnological fit to the data. More precisely 
there is an initial exponential growth (until 1-3 mm in di- 
ameter) followed by the vascular Gompertzian phaseQ- 
Then it seems reasonable to think that cancer growth fol- 
lows a general pattern that one can hope to describe by 
macroscopic variables, constrained by a set of environ- 
mental conditions crucial to understand the fundamental 
features of the growth as, for example, the shape and the 
maximum possible size of the tumor and the onset of the 
tissue invasion and metastasis. Following this line of re- 
search, for example, the universal model proposed in Q 
has been recently applied to cancer^. 

In this paper we consider a macroscopic model of tu- 
mor growth that: i) gives an energetic basis to the Gom- 
pertzian law; ii) clearly distinguishes among the gen- 
eral evolution patterns, which include feedback effects 
and external constraints; iii) can give indications on the 
different tumor phases during its evolution. The pro- 
posed macroscopic approach is not in competition with 
microscopic models |5j, but it is a complementary in- 
strument for the description of the tumor growth. The 
Gompertzian curve is solution of the equation 



dt ' V N 
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where N(t) is the cell number at time t, 7 is a constant 
and is the theoretical saturation value for t — > 00. 

It is quite natural to identify the right hand side of 
Eq. 0} as the number of proliferating cells at time t and 
then to consider f p (N) = 7 In (-#*■) as the fraction of 
proliferating cells and 1 — f p (N) — f np the fraction of 
non proliferating cells. Here, we observe the difference 
with respect to an exponential growth that corresponds 
to a N independent specific proliferation rate, typical of 



a system of independent cells. Instead, the logarithmic 
dependence observed in the Gompertzian law, can be re- 
garded as a sort of feedback mechanism, such that the 
system, at any time, rearranges the growth rate and the 
fraction of proliferating cells according to the logarithm 
of the total number of cells N(i). 

Living tumor cells consume an amount of energy, un- 
der the form of nutrient and oxygen, which is recovered 
from the environment (which will be addressed in the fol- 
lowing as B), represented by a living body (in vivo) or 
an external solution (in vitro). B, as already mentioned, 
has also a mechanical action that can be responsible of 
compression, increase of pressure and density and change 
of shape of the tumor (addressed in the following as A). 
Starting from these observations, we intend to describe 
the tumor growth on an energetic basis. 

Let us first consider the problem at fixed time and re- 
call the first law of thermodynamics. A change of the 
energy content of a system has three possible origins, 
namely mechanical work performed on or by the system, 
variation of the number of elementary constituents (par- 
ticles), heat exchange with the environment. Accord- 
ingly, the internal energy of the system is expressed as the 
sum of three terms, respectively —PV (minus the prod- 
uct of pressure times volume), fj,N (where fj,, the chemi- 
cal potential, is the change in energy when the number of 
particles N is increased by one unit), S/(3 (heat content). 

The basic ingredient in our model is an energy bal- 
ance law for the tumor in the same spirit of the first 
law of thermodynamics. However, it should be clear that 
the variables that we use to describe our system are not 
properly those introduced in thermodynamics but rather 
the rules and the relations among these variables are the 
same as those of thermodynamics. In fact at some fixed 
time t one can regard the available energy for the tumor, 
U, as the sum of three pieces, U = fl + fiN + Em, where 
il = —PV is associated to the mechanical work necessary 
to contrast the external pressure and to cause a change 
of shape, /iN is proportional to the number of cells N of 
A and fi represents the variation in energy related to a 
change in N (note, for the sake of clarity, that the me- 
chanical work related to the change in volume of A due 
to the increasing number of cells is taken into account in 
f2) and finally Em , which represents the energy spent for 
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metabolic activities. Since the change in the number of 
cells in the tumor is due to biological proliferation and 
by recalling the above definition of fiN, it is natural to 
associate fj,N to the energy that A needs to generate one 
new cell, i. e. the energy necessary for reproduction, that 
has to be distinguished from the energy associated to all 
the other metabolic processes which is addressed as Em ■ 

A set of macroscopic variables suitable to describe A 
is: the volume of the whole tumor, V, the variable /i 
introduced above that we address as chemical potential, 
in analogy with the thermodynamical quantity, and (3, 
analogous to the thermodynamical inverse temperature, 
which is introduced, together with its conjugate variable 
S (that corresponds to the thermodynamical entropy), 
by the relation Em = S/f3. The variables V, [i and (3 
define the macroscopic state of A at time t. 

Statistical mechanics provides a framework that ex- 
plains the basic laws of thermodynamics starting from 
a microscopic description of the thermodynamical sys- 
tem. Therefore, having in mind the energetic aspect of 
the problem, we shall introduce a minimal microscopic 
model of the cell, where the biological aspects are re- 
duced to the indispensable, which however catches the 
essential ingredients that provide a sound macroscopic 
description and a good comparison with the quantitative 
data. In an extremely simplified picture, we assume that 
the microscopic states of each cell can be labelled by an 
integer index I starting from a ground state level with 
energy e so that each cell has an energy spectrum of the 
form ti = e + Id, where 5 is the minimum positive en- 
ergy gap between two states. It should be clear that here 
the energy of a cell indicates the energy that the cell is 
consuming per unit time, so that it is obvious that the 
minimum energy in the spectrum must be positive e > 0. 

According to the basic rules of statistical mechanics, 
once we have the energy spectrum of a single cell we 
can determine the grand partition function, Z((3, V, /i) = 
II^ exp (e~^^ l ~^) and the grand potential, correspond- 
ing to the quantity Q introduced above: 

mV,n) = -^Z = -je-^-^ (2) 

where W = 1/(1 — e _/3<5 ). From the derivatives of 
V,n) we respectively get N, for constant V and (3, 
and Em, for constant V and /x: 




where C = W (38 exp(-/3(5). From Eqs. © and © 
it is straightforward to express /i in terms of N : [i = 
e+ (1/(3) \n(N/W) and also to derive the average energy 



per cell: U/N = (E M + (iN + Q)/N = e + C/(3, which 
shows that U/N depends on one macroscopic variable 
only, the inverse temperature (3, as it is the case for an 
ideal gas in thermodynamics. 

At this point we observe that there are two natural 
constraints in our model: in fact it is clear that an in- 
crease in the number of cells N requires a positive energy 
supply which has a direct consequence on the sign of the 
chemical potential: \x > 0. Moreover Em > because 
the cells also require a positive energy supply for their 
metabolic activity. In addition, in order to compensate 
any mechanical stress on the tumor we need to require 
a stronger constraint: Em + ^ > 0. According to the 
relations derived above, the constraint /i > implies a 
minimum number of cells TV < N m = W exp(— (3e) and 
(E M + n)/N = (1/(3) (C + \n(W/N)) > provides a 
maximum: N < Nrx, — N m exp(C + (3e). Then, by re- 
placing N m and in the previous relations we get fj, = 
(l/P)\n(N/N m ) and (E M + il)/N = (1/(3) ln^/jV). 
Provided one identifies the maximum number of cells here 
obtained with the asymptotic value of N appearing in Eq. 
(JIJ, the latter equation for (Em + Cl)/N is remarkably 
similar to the the fraction of proliferating cells at time 
t, f p (t), introduced after Eq. © and this will be our 
starting point to derive the Gompertzian law. 

Now we assume that A, due to cellular biological du- 
plication and possible changes induced by B, can modify 
its macroscopic state and therefore also the macroscopic 
variables which define the state itself and determine the 
grand partition function Z. We assume that these changes 
occur on a time scale that is much larger than the typi- 
cal time during which the macroscopic variables reach the 
equilibrium and become constant, so that we can think of 
the evolution of A as a slow change from one macroscopic 
state, where the statistical framework is well defined, to 
another one, where the same framework, but with differ- 
ent values of the macroscopic variables, holds. 

Recalling that at any time fixed t one can write the 
identity 1 = (f p + f np ) and the energy sum rule is 
U = (C/(3 + e)(f p + f np )N = E M + Q + fiN and, ac- 
cording to the observed similarity between the expres- 
sions of f p and (Em + ty/N, it looks natural to iden- 
tify E M + ft — (C I ' (3 + e)f p N and, correspondingly 
fi = (C//3+ e)f np This basic assumption is supported by 
a biological counterpart. In fact in |6( the feedback effect 
of Eq. © is described by introducing some biochemical 
effects due to growth inhibitor factors, released by the 
dead cells, which inhibit the proliferation and eventually 
are responsible of the saturation of the tumor size, so 
that the concentration of these inhibitor factors should 
be proportional to f np . In our energetic picture the con- 
centration of inhibitor factors should act enhancing the 
energy necessary for increasing, through cellular dupli- 
cation, the number of cells in the tumor and therefore 
enhancing /i. Then, the presence of growth inhibitor fac- 
tors leads to a qualitative proportionality between f np 
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and \i that supports our assumption. 

With the statement E M + ?L = (C//3 + e)f p N it is 
straightforward to derive the Gompertz equation. In 
fact, by definition, AN in an interval At is proportional 
to the number of proliferating cells and then AN = 
c\At f p (t)N — c^Af (Em + ^) with c%, C2 constants. 
When At — > one recovers Eq. |JTJ with 7 = c 2 //3. Re- 
markably, this model explains the linear correlation be- 
tween the two standard parameters used to fit the data 
on tumor growth by the Gompertzian curve 0, 0] ■ 

A comment is in order. N m is related to ji = 0. In 
our set of relations N < N m corresponds to fi < which 
is not physically allowed. This indicates that our model 
and the feedback mechanism embedded in it, start being 
valid and hold only when N > N m . This again fits with 
the experimental observations that the onset of the Gom- 
pertzian growth does not occur at the beginning of the 
tumor development but only after an exponential phase 
growth . It is reasonable to think that when N < N m , 
due to the abundance of energy provided by B, the cells of 
A duplicate as a system of independent cells with no fur- 
ther energetic effort for the duplication due to the pres- 
ence of the other cells, i.e. with \i = (and not \i < 0), 
which indicates that the feedback effect should not take 
place and, as noticed after Eq. this implies an ex- 
ponential growth. Therefore, even if our model does not 
work for N < N m , one reasonably expects a continuous 
crossover between exponential and Gompertzian growth 
at TV = N m . 

To analyze some phenomenological implications of our 
model we shall consider in the following the simple case 
of the multicellular tumor spheroids (MTS), which follow 
the Gompertzian growth law as tumors in vivo, despite 
they have an avascular growth. We shall first consider 
experiments without external stress and, later, the effects 
due to a change of the external pressure. 

i) ENERGETIC GROWTH 
Here we shall refer to experiments where the variation of 
the external supply of nutrient and/or oxygen is studied 
0, Q and we shall rely on a minimal description of the 
MTS where |2j: a) the thickness, k, of the layer where 
the nutrient and oxygen is delivered (the crust) is inde- 
pendent on the spheroid radius R (k as a function of R is 
constant within 10%); b) the cellular density is constant; 
and essentially independent of the nutrient concentra- 
tion; c) the difference in nutrient and oxygen supply for 
the cells in the inner and outer layers of the crust is neg- 
ligible in first approximation, d) the cells at distances 
d < R — k ii R > k from the center of the spheroid are 
dead and only the cells of the crust are alive. 

The data, for various concentrations of glucose and 
oxygen show that 6]: 1) after an initial exponential 
growth (3-4 days) where the doubling time is indepen- 
dent of the oxygen and glucose concentration, the radius 
R of MTS grows according to a Gompertzian law as the 
one in Eq. Q , where N is replaced by R and Noo by Roo 
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FIG. 1: Thickness (fim) vs. glucose concentration (mM). 
Figure (a) is for an oxygen concentration of 0.28 mM and 
obtained for AT°, = 2.1 10 5 , a = 1.8 10 5 , G° c = 0.8. Figure (b) 
is for an oxygen concentration of 0.07 mM and obtained for 
iV^ = 7.6 10 3 , a = 1.03 10 4 , G° = 0.8. 

(Roo is the maximum radius that the MTS reaches at sat- 
uration at very large i); 2) there is an approximate linear 
correlation between k and the logarithm of the glucose 
concentration, G c , at fixed oxygen quantity; 3) there is 
a correlation between the radius at the onset of necrosis, 
R = k and at the growth saturation, R = Roo. 

It is reasonable to think that at the early stages of the 
growth, when R < k, all cells receive the nutrient and 
oxygen supply while, when R > k, due to a more diffi- 
cult diffusion of the nutrient, there is a fraction of non- 
proliferating cells and the feedback effect starts. Accord- 
ing to our model, one observes the Gompertzian phase 
for t > t*, with t* defined by the condition R(t*) = k, 
when the feedback begins. Since there is a simple ge- 
ometrical relation between N and R due to point b), 
it is easy to check that, for R >> k, our model pre- 
dicts a Gompertzian growth for the MTS radius up to 
an asymptotic value Roo, which is determined by Noo- 
In our model, the number of cells at time t* is fixed by 
the condition fj, = 0, i.e. N(t*) = N m and, as seen be- 
fore, N m — A oc exp(— e/3 — C). Putting all together these 

1 /3 

informations, we find k = R(t*) ~ No£ . 

Since the experimental data at fixed oxygen quantity, 
are almost consistent with a linear relation Noo = a(G c — 
G°) + where a, G° and N° are ad hoc constants (see 
Figure captions for details), jal^' we are arj le to relate 
the MTS radius at the onset on necrosis, k, to G c : 

k(G c ) = a (NU 3 (G C ) - N°J /3 ) + k (5) 

where a and ko are constants depending on the supplied 
oxygen. From Eq. (J5J one obtains the correlation among 
Noo, G c and k. In Fig.l and Fig. 2 the previous behav- 
iors are compared with data without optimization of the 
parameters. Despite the oversimplified approach, there 
is a good agreement with the data. Moreover the pro- 
posed model predicts the Gompertzian behavior on the 
basis of a feedback effect and a constant supply of nutri- 
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FIG. 2: Spheroid saturation cells number vs. diameter (/oti) 
at which necrosis first develops. Circles refer to culture in 
0.28 mM of oxygen and the relative curve is obtained for 
a — 1.78 and 2ko = 216 [im. Triangles refer to culture in 
0.07 mM of oxygen and in this case a — 5.3, 2fco = 100 fim. 



ent for cell and this explains why, under these conditions, 
the MTS and the in situ tumors follow the same general 
growth law. 

ii) BIOMECHANICAL EFFECTS 
The heterogeneous elastic characteristics of the environ- 
ment give the strongest external constraints[9j. Again 
the MTS are useful to study how the mechanical fea- 
tures of the environment can modify the growth pattern. 
When the MTS are under a solid stress, obtained for ex- 
ample by a gel, the experimental data indicate that 

1) an increase of the external gel concentration C g in- 
hibits the growth of MTS reducing the size of Roo with 
respect to the case with C g = (as before is the 
maximum radius that the MTS reaches at saturation); 

2) the cellular density at saturation increases with 

In our model the mechanical energy is included in the 
energetic balance by the term f2, and therefore in this 
problem f2 depends on C g . In particular, for growths 
in absence of gel, C g = 0, we neglect its value, ft ~ 0, 
because of the very low pressure P ~ 0, if compared to 



the case where C g ^ 0. Also Noo depends on Q and, 
by recalling that Aoo(f2) is obtained from the constraint 
Em + > and therefore when f2 = one can get 
Noo(0) directly from Em > 0, it is straightforward to 
derive A^fi) = Aoo(0)/e. For a growth with constant 
density this reduction should also imply a decrease of the 
maximum size of the spheroids, i.e. i?oo(r2) < i?oo(0). 
For a quantitative approach one has to take into account 
the correlation between the density p and the dynamical 
growth and, since SI = —N/j3 = —PV, the pressure is 
P(t) = p(t)/(3. So, for poo, one has [poc(ty - Poo(0)] = 
(3P e where P e is the external pressure due to the gel C g , 
and poo(0) is the density for P e = 0, i.e. for C g = 0. 
Then increases with C g according to point 2). A lin- 
ear fit to the data gives Poo(Jl) — Poo(0) = 0.71C g @, and 
since p ~ N/R 3 , one still needs A^f^/A^O) to obtain 
the dependence of the MTS radius at saturation i? OQ (f2) 
on C g . The relation N^Cl) — iVoo(0)/e, derived above 



is discontinuous in the limit r2 — > 0, which is an artifact 
of our oversimplified microscopic model. However, for 
intermediate values of the gel concentration we do not 
expect large deviations from this result and we use it, 
in first approximation, to determine i?oo(f2) in terms of 
C g . According to this procedure we have no adjustable 
parameters left and the comparison with the data is re- 
ported in Table I for C g in the range 0.3 — 0.8 %. 



C g (percent) 



0.3 
0.5 
0.7 
0.8 



2_Roo(Sl) [fJ.m] exper. 



450 
414 
370 
363 



2R oo {0) [fim] fit 



452 
429 
404 
394 



TABLE I: Comparison with the experimental data as dis- 
cussed in the text. The experimental error is about ±10% 
(see Fig. la in ||). 



In conclusion, in the presented model the Gompertzian 
growth is essentially due to an average energetic balance 
which represents a macroscopic description of a complex 
underlying dynamics. The comparison with the MTS 
data suggests that it is a good starting point for a macro- 
scopic analysis of tumors in vivo and indeed the Gom- 
pertz equation has interesting properties that can explain 
why this approximation works in realistic conditions |l0l| . 
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